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ABSTRACT 

We describe the radio evolution of SN 2001ig in NGC 7424, from 700 days of 
multi-frequency monitoring with the Australia Telescope Compact Array (ATCA) 
and the Very Large Array (VLA). We find that deviations of the radio light curves 
at each frequency from the standard "minishcll" model are consistent with density 
modulations in the circumstellar medium (CSM), which seem to recur with a period 
near 150 days. One possibility is that these are due to enhanced mass-loss from thermal 
pulses in an AGB star progenitor. A more likely scenario however is that the progenitor 
was a Wolf-Rayet star, whose stellar wind collided with that from a massive hot 
companion on an eccentric 100 day orbit, leading to a regular build-up of CSM material 
on the required time and spatial scales. Recent observations of "dusty pinwheels" in 
Wolf-Rayet binary systems lend credibility to this model. Since such binary systems are 
also thought to provide the necessary conditions for envelope-stripping which would 
cause the Wolf-Rayet star to appear as a Type Ib/c supernova event rather than 
a Type II, these radio observations of SN 2001ig may provide the key to linking 
Type Ib/c SNe to Type lib events, and even to some types of Gamma-Ray Bursts. 

Key words: circumstellar matter - supernovae: individual: SN 2001ig - binaries - 
stars: Wolf-Rayet - galaxies: individual: NGC 7424 - gamma rays: bursts. 



1 INTRODUCTION 

Radio studies of supernovae (SNe) can provide valuable in- 
formation about the density structure of the circumstellar 
medium, the late stages of stellar mass-loss, and indepen- 
dent distance estimates for the host galaxies dWeiler et alJ 
2002). Furthermore, with the growing realisation that 
some Gamma-Ray Bursts may be intimat ely linked with 
SNe [e.g., GRB980425 and SN 1 998bw jKulkarni et alJ 
Il998lh GRB011121 and SN 2001ke iGarnavich et al.ll2003lh 
GRB030329 and SN 2003dh JStanek et all2003l) 1. such stud- 
ies are crucial to understanding GRB environments. To 
date, radio emission has only ever been detected from core- 
collapse SNe of Type II and Type Ib/c [for a recent review 
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of supernova taxonomy, see iTurattd i200 3)l . and not at all 
from thermonuclear Type la SNe. 

Supernova 2001ig was discovered visually by lEvansI 

J200ll) on 2001 Dec 10.43 UT in the outskirts of the nearby 
late-ty pe spiral galaxy NGC 7424 (D = 11.5 Mpc; iTullvl 
( 1988)). No supernovae have been recorded previously in this 
galaxy. Initial optical spe ctroscopy of SN 2 001ig from Las 
Campanas Observatory iMatheson fc Jh al l200_j) reveale d 
similarities to the Type lib SN 1987K l|FiliDPenkolll988h . 
while spectra from th e European Southern Observatory over 
the fo llowing month llClocchiatti fc Prietdl200ll : IClocchiatti 
2002) showed a similar behaviour to that of the Type lib 
SN 1993J, by transitioning from Type II to Type Ib/c as the 
H recombination lines weakened. By October 2002, the tran- 
sition to a Type Ib/c in the nebular pha se was well and truly 
complete ( Fili ppenko fc Chornockl20'02Tl . SN 2001ig was also 
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detected by the ACIS-S instrument on board the Chandra 
X-ray Observatory on 2002 May 22 UT with a 0.2- 10.0 keV 
luminosity ~ 10 38 erg s" 1 jSchleeel fc Rvdeill2002Tl . A sec- 
ond observation three weeks later showed that t he X-ray 
luminosity had halved in that time dSchlegeJl20oil . 

Radio monitoring of SN 2001ig with the Australia Tele- 
scope Compact Array 1 (ATCA) commenced within a week 
of its discovery, and has continued on a regular basis. In Sec- 
tion |3| we present multi-frequency radio flux data from the 
first 700 days. We describe our attempts to fit the radio 'light 
curves' with a circumstellar interaction model in Section [I] 
and discuss the deviations from this model in more detail in 
Section |5] Our conclusions are presented in Section |5] 



2 RADIO MONITORING 

Table^contains the complete log of radio flux measurements 
from the ATCA. Column 2 lists the days elapsed since ex- 
plosion, which is derived from the model fitting in Section 0] 
Total time on-source ranged from as little as 2 hours, up to a 
full 12 hour synthesis, but was typically 4-6 hours. Since the 
ATCA is capable of observing in 2 frequency bands simulta- 
neously, determining fluxes in 4 frequency bands on the same 
day required time-sharing. Dual-frequency observations cen- 
tered on 18.75 and 18.88 GHz and bandwidths of 128 MHz 
were carried out using a prototype receiver system on just 
three ATCA antennas. The central frequencies of the other 
bands are 8.640, 4.790, 2.496, and 1.376 GHz, and the band- 
width is 128 MHz. From 2002 July 11 onwards, the S-band 
central frequency was changed from 2.496 to 2.368 GHz, to 
reduce the amount of in-band interference. The ATCA pri- 
mary flux calibrator, PKS B1934-638 has been observed once 
per run, while observations of the nearby source PKS B2310- 
417 allow us to monitor and correct for variations in gain and 
phase during each run. 

The data for each observation and frequency have been 
edited and calibrated using the miriad software package. 
Rather than image the visibility datasets, then clean them 
to some arbitrary level, the UVFIT task is used instead in 
the visibility domain to fit simultaneously a point source at 
the known location of SN 2001ig, as well as a background 
source fortuitously located just 18.5 arcsec to the south- 
west (Figure EJ. This source, located at a = 22 h 57 m 29.6 s , 
S = -41°02'40" (J2000) does not appear in any radio 
source catalogue, but was found to have the following 
fluxes as at 2001 Dec 15 UT: ,5(8.640 GHz) = 5.8 mjy, 
5(4.790 GHz) = 12.9 mjy, 5(2.496 GHz) = 26.4 mjy, and 
5(1.376 GHz) = 47.0 mjy. After scaling the background 
source to these values at each epoch, the SN 2001ig fluxes 
have been scaled accordingly. The uncertainties in Table 
include both the formal fitting errors, and the possibility 
that this background source is intrinsically variable (as is 
quite likely at the higher frequencies) . This technique of fit- 
ting the SN flux density in the visibility domain, then boot- 
strapping to the flux of the adjacent source, was crucial to 
recovering valid flux measurements from observations with 
poor phase stability and/or limited hour-angle coverage. 

The Australia Telescope is funded by the Commonwealth of 
Australia for operation as a National Facility managed by CSIRO. 



The ATCA observations have been supplemented by a 
few early observations with the Very Large Array (VLA) 2 . 
The obse rving and data anal ysis procedure follows that de- 
scribed in Wei ler et al] Il986t> , and the results are listed sep- 
arately in Tabled Owing to the low elevation of SN 2001ig 
as observed from the VLA, and the compact configuration 
of the VLA at the time, both the sensitivity and the resolu- 
tion are relatively poor. Nevertheless, these data prove to be 
important in constraining the early evolution of SN 2001ig, 
particularly at the very highest frequencies. 



3 RADIO LIGHT CURVES 

The ATCA and VLA radio data is plotted in Figure The 
data at 15.0 GHz are not shown here, to reduce confusion 
with the other high-frequency points, but is incorporated 
in the model fitting of Section 2] The time evolution of the 
spectral index a (where flux S oc v +a ) between simultaneous 
positive detections at 1.4 and 4.8/4.9 GHz, and between 
4.8/4.9 and 8.6 GHz, is plotted in Figure El 

The radio "light curve" of a Type II supernova can be 
broadly divided into three phases: first, there is a rapid turn- 
on with a steep spectral index (a > 2, so the SN is brightest 
at the higher frequencies), due to a decrease in the line-of- 
sight absorption. After some weeks or months have elapsed, 
the flux reaches a peak, turning over first at the highest fre- 
quencies. Eventually, the SN begins to fade steadily, and at 
the same rate at all frequencies, in the optically-thin phase. 

Although broadly consistent with this picture, the radio 
light curve of SN 2001ig displays significant departures from 
a smooth turnover and decline, which are most pronounced 
at 8.64 and 4.79 GHz. At around day 80, the flux at these 
frequencies reversed its initial decline, and by day 130 had 
almost doubled. The flux remained almost constant for a pe- 
riod of ~ 50 days, before resuming its decline at close to the 
original rate. Near day 250, the decline was again temporar- 
ily interrupted for another 50 days. There are indications of 
perhaps one more bump after day 450, but by this stage the 
SN has faded to the milli-Jansky level where any variations 
are comparable to the measurement uncertainties. The de- 
viations are less pronounced, but still evident at 2.50 and 
1.38 GHz. 



4 MODEL FITS 

The general properties of supernova radio light curves as 
outlined above are quite well re presented by a mo dified ver- 
sion of the "minishell" model of IChevalierl £l982), and have 
been successfully parameterise d for m ore than a dozen RSNe 
(see Table 2 of IWeiler et all fc!002h L Radio synchrotron 
emission is produced when the SN shock wave ploughs into 
an unusually de nse circumste ll ar me dium (CSM). Following 
the no tation of IWeiler et al] (I2002T) and ISramek fc Weilerl 
(2003), we model the multi- frequency evolution as: 



2 The VLA telescope of the National Radio Astronomy Observa- 
tory is operated by Associated Universities, Inc. under a cooper- 
ative agreement with the National Science Foundation 
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Table 1. SN 2001ig radio flux measurements with the ATCA. 



Date 


Days since 


5(18.8 GHz) 


5(8.6 GHz) 


5(4.8 GHz) 


5(2.4 GHz) 


5(1.4 GHz) 


(UT) 


2001 Dec 3 UT 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


15/12/01 


12 




2.1 ± 0.3 


0.6 ± 0.3 




< 0.8 


18/12/01 


15 




3.9 ± 1.6 


1.6 ± 0.2 






99 /1 9 /01 


19 




8.7 ± 0.8 


2.7 ± 0.3 






26/12/01 


23 




15 + 30 


4.8 ± 0.3 






31 /12/01 


28 


43 ± 4 




6.3 ± 0.3 


< 4.6 


< 1.8 


07/01 /09 


35 






10 6-1-03 
J-U . u in u.o 






1 /01 /09 


38 


22 ± 3 


23.6 ±4.1 


14.2 ± 0.9 






1 ^ /01 /09 


43 


11 ± 4 


ro . y zc o . u 


15.1 ± 0.7 






90/01 /09 


48 








5.3 ± 0.3 




09 /09 /09 

UZ / UZ J OA 


61 




9.9 ±22 


18 4 + 44 






1 7/09 /09 


76 




7.0 ± 1,1 


11.0 ± 0.5 






9fi /09 /09 

ZU / UZ / UZ 


85 








11.5 ± 0.3 


7.0 ± 0.5 


1 7 /0^/09 


104 




10.7 ± 2.0 


18.1 ± 1.3 






1 9 /03 /02 


106 








25.1 ± 1.3 


13.8 ±1.5 


28/03/02 


115 




11.1 ± 2.0 


21.9 ± 1.0 


23 5 + 3 


14.2 ± 0.6 


09/04/02 


127 




12.6 ± 0.7 


oi O + 06 






99 /OzL/09 


140 




12.7 ± 0.9 


21.6 ± 0.7 






12/05/02 


160 




12.4 ± 0.5 


21.7 ± 0.4 


30 3 + 1 6 


21.8 ± 1.0 


27/05/02 


175 




8.0 ± 1.7 


15.6 ± 1.3 


24.4 ± 0.3 


21.5 ± 0.7 


09/06/02 


188 




6.2 ± 1.7 


11.9 ± 1.2 


19.7 ± 1.0 


19 5 + 06 


30/06/02 


209 




4.6 ±1.4 


8.9 ± 1.0 


14.4 ± 0.6 


in k _|_ o a 

J- CI • <J _1_ \J • <J 


1 1 /07/02 


220 




4.0 ± 1.2 


7.6 ± 0.9 


14.9 ± 0.4 


19.8 ± 1.2 


28/07/02 


237 


- 


3.6 ±0.9 


7.1 ± 0.7 


13.0 ±0.3 


18.2 ±0.5 


11/08/02 


251 




3.3 ± 1.0 


6.4 ± 0.8 


12.7 ±0.3 


18.1 ± 1.0 


23/08/02 


263 




3.8 ±0.5 


6.7 ± 0.3 


12.5 ±0.2 


17.8 ±0.9 


31/08/02 


271 




3.7 ±0.3 


6.7 ± 0.2 


12.5 ±0.2 


17.4 ±0.9 


17/09/02 


288 




3.2 ± 1.0 


6.6 ±0.7 


12.3 ±0.3 


16.4 ±0. 6 a 


12/10/02 


313 




2.7 ±0.9 


6.7 ± 0.4 


10.2 ±0.6 


15.4 ±0.8 


22/11/02 


354 




2.5 ± 0.3 


4.5 ±0.3 


8.0 ±0.3 


11.7 ±0.5 


17/12/02 


379 




1.3 ±0.3 


3.6 ±0.3 


6.5 ±0.5 


11.5 ±1.1 


05/02/03 


129 




1.8 ±0.4 


2.2 ±0.2 


6.1 ±0.5 


10.0 ±1.3 


16/03/03 


468 




1.5 ± 0.4 


2.6 ±0.3 


5.5 ±0.4 


8.8 ±0.4 


20/05/03 


533 




1.2 ± 0.5 


2.4 ± 0.2 


4.5 ±0.2 


7.7 ±0.7 


03/08/03 


608 




0.9 ±0.4 


2.0 ±0.2 


3.8 ±0.3 


6.3 ±0.4 


07/11/03 


704 




0.8 ±0.3 


1.6 ±0.2 


3.1 ±0.3 


5.1 ±0.5 



a IChandra fc Ravi i2002ft reported a flux at 1.4 GHz of 11.7 ± 1.5 mjy for SN 2001ig on 2002 Sep 25.8 UT 
with the Giant Meterwave Radio Telescope (GMRT). T he source of this discre pancy has since been traced 
to an elevation-dependent gain error in the GMRT data ( Chandra & Rav 2003 ) . 



Table 2. SN 2001ig radio flux measurements with the VLA. 



Date 


Days since 


5(22.5 GHz) 


5(15.0 GHz) 


5(8.5 GHz) 


5(4.9 GHz) 


5(1.4 GHz) 


(UT) 


2001 Dec 3 UT 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


28/12/01 


25 


37.2 ±5.6 


36.5 ±3.7 








07/01/02 


35 


22.3 ±3.4 


30.9 ±3.1 








10/01/02 


38 


16.6 ± 2.5 


25.1 ± 2.6 








13/01/02 


11 


15.8 ± 2.4 


20.5 ±2.1 


23.5 ± 1.8 






17/01/02 


45 










1.6 ±0.4 


27/01/02 


55 








9.8 ± 0.7 


3.2 ±0.3 


21/03/02 


108 






10.5 ±0.5 


18.6 ± 0.9 


12.6 ±0.8 



with 



Toxtcrnal = TCSM homog + Tdistant (2) 

where 
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Figure 1. Contours of 4.790 GHz radio emission, on a V-band image of NGC 7424. The radio observations were made with the ATCA on 
2002 Feb 17 UT, and have a synthesised beamwidth of 3.8 X 1.5 arcsec. SN 2001ig is the upper-left of the two (unres olved) sources. The 
optic al image is 3.5 arcmin on a side, obtained with the DFOSC instrument on the Danish 1.54-m telescope on La Silla llLarsen fe Richtlerl 
1999), and made available through NED. 



TCSM homog =^ 2 (^) (| 

Tdistant = Mri^) ' 

and 



to 



day 



rcsM clumpo =if3(rik) (I 



t -t 
day 



(3) 
(4) 

(5) 



with the various K terms representing the flux density (Ki), 
the attenuation by a homogeneous absorbing medium (K2, 
K4), and by a clumpy/filamentary medium (K3), at a fre- 
quency of 5 GHz one day after the explosion date to- The 
TcsM homog and 7~csM chlmp3 absorption arises in the circum- 
stellar medium external to the blast wave, while Tdistant is a 
time- independent absorption produced by e.g., a foreground 



Hn region or more distant parts of the CSM unaffected by 
the shock wave. The spectral index is a; j3 gives the rate 
of decline in the optically-thin phase; and 5 and 5' describe 
the time dependence of the optical depths in the local ho- 
mogeneous, and clu mpy / filamentary CSM ; respec tively (see 
IWeiler et al.1 <2002t and ISramek fc Weilerl <2003ft for a de- 
tailed account of how these parameters are related). For 
lack of sufficient high-frequency data prior to the turnover 
to constrain it, we adopt internal = 0. 

In order to assess the gross properties of SN 2001ig in 
comparison with other Type lib RSNe, we have fit this stan- 
dard model to the data in Tables Q and D2J but excluding 
days 48-70, and days 110-190. Thus, the model fit is con- 
strained primarily by the rise at early times, the region of 
the high-frequency turnover, and by the late-time decay. The 
actual date of explosion to is found to be 2001 Dec 3 UT, 
one week prior to discovery. The full set of model parameters 
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T 1 1 f 




(t-t )/days [t = (10-Dec-01) -7.194 days.] 

Figure 2. SN 2001ig at radio frequencies of 22.5/18.8 GHz (circles, thick solid line), 8.6/8.5 GHz (crosses, dashed line), 4.9/4.8 GHz 
(squares, dash-dotted line), 2.4 GHz (triangles, dotted line), and 1.4 GHz (diamonds, dash-triple dotted line). The curves are a model fit 
to the data, as described in the text. 



which yields the minimum reduced \ 2 is given in Table [31 
and the model curves are plotted in Fig. [21 For comparison, 
we show in Table [5] the equivalent paramete rs for two other 
well-s ampled Type lib RSNe: S N 1993J llVan Dvk et"afl 
12004) in M81, and SN 2001gd jStockdale et all l2003h in 
NGC 5033. Note tha t we have fixed the value of S to be 
(a — (3 — 3), as in the IChevalierl il982f) model for expansion 
into a CSM with density decreasing as r~ 2 . 



The spectral index a of SN 2001ig is virtually identi- 
cal to SN 1993J, but less steep than SN 2001gd. However, 
the rate of decline f3 is much steeper in SN 2001ig than in 
either of the other Type lib's, and the time to reach peak 
5 GHz flux is also much shorter. In this respect, SN 2001ig 
has behaved more like a Type Ib/c SN than most "nor- 
mal" Type II SNe. The interpolated peak 5 GHz luminosity 
would be about twice that attained by SN 1993J, though in 
practice SN 2001ig was near a local minimum in the flux at 
that time, and the actual peak was not reached for another 
40 days. 



Using the methodology outlined in lWeiler et alJ J2002t) 
and ISramek fc Weilerl i2003l) . we can derive an estimate of 
the progenitor's mass-loss rate, based on its radio absorp- 
tion properties. Substituting our model fit results above 



into their equation 11, and assuming both TcsM homog and 
7"csM clumps contribute to the absorption 3 , we find that 

^ =(2.2 ±0.5) x 10- 
w/(10 km s" 1 ) v ; 

where w is the mass-loss wind velocity, and the ejecta 
velocity as measured from the earliest o ptical spectra 
JClocchiatti fc Prietdl200lHClocchiattll2002l) is in the range 
15,000 - 20,000 km s _1 . Clearly, this is only an average value, 
subject to major variations discussed in the next section, but 
is in the same domain as the mass-loss rates derived similarly 
for SN 1993J and SN 2001gd (Table 0. Though generally 
less well-constrained, mass-loss rates may also be estimated 
directly from the radio emission properties, relying only on 
the peak 5 GHz luminosity, and the time tak e n to reach 
that peak. Equations 17 and 18 of IWeiler et alJ i2002l) give 
M/w = 1.5 x 10" 5 and 3.5 x 10~ 5 for the average Type Ib/c 



3 Case 2 in the notation of IWeiler et all <2002l) . Note that 
there is a misprint in that section, namely that in the limit of 
r CSM clumps -> 0, then (t°£) -» T"cSM homoK in equation 13, and 
not r CSM homog - We prefer the term "homogeneous" here over 
"uniform" , since the latter could give the misleading impression of 
no density gradient at all, whereas an r~ 2 dependence of density 
is implicit. 



6 S. D. Ryder et al. 



02 



X 
CD 



CO 

■+J 
o 

CD 
CO 



c\2 

I 



T 



T 



• 1.38 - 4.79 GHz 
□ 4.79 - 8.64 GHz 



H ♦- 



+ 

H[] [][] [] [] [][] 



o 



200 



400 



600 



time [days since explosion] 



Figure 3. Evolution of spectral index a for SN 2001ig, plotted linearly as a function of time, between 1.38 and 4.8/4.9 GHz (solid 
circles); and between 4.8/4.9 and 8.64 GHz (open squares). 



Table 3. Comparison of radio light curve model parameters. 



Parameter 




SN 2001ig 


SN 1993J 


SN 2001gd 






2.71 x 10 4 


1.36 x 10 4 


1.49 x 10 3 


a. 




-1.06 


-1.05 


-1.38 


P 




-1.50 


-0.88 


-0.96 


Kz 




1.38 x 10 3 


9.14 x 10 2 


3.25 x 10 6 


5 




-2.56 


-1.88 




K 3 




1.47 x 10 5 


8.33 x 10 4 


1.05 x 10 3 


5' 




-2.69 


-2.26 


-1.27 


K A 




0.0 


2.76 x 10" 3 




Time to L 5 GHz peak 


(days) 


74 


167 


173 


Lb GHz peak (erg s" 1 


Hz- 1 ) 


3.5 x 10 27 


1.4 x 10 27 


2.9 x 10 27 


Mass-loss rate (Mq 


yr" 1 ) 


(2.2 ± 0.5) x 10~ 5 


2.1 x 10" 5 


3.0 x 10" 5 



SN and Type II SN, respectively. Thus, the mass-loss rate 
calculations are in good agreement, with SN 2001ig interme- 
diate between the expected rates for Type Ib/c and Type II 
SNe, consistent with its Type lib classification. 



5 DISCUSSION 

In Figure 2] we have plotted the deviations of the observed 
flux density, from the best-fit model curves as shown in Fig- 
ure |21 The solid line in this figure is a 4-point boxcar aver- 
age of the mean deviation over all frequencies at each epoch, 
which serves to emphasise the quasi-damped harmonic na- 
ture of the deviations, having a period near 150 days, and 
peak intensity declining with time (i.e., with increasing dis- 
tance from the star). The r.m.s. deviation of the actual data 



from the smoothed interpolation is less than 1/3 of the am- 
plitude of the observed modulation. Since the fractional am- 
plitude of these deviations is virtually identical at each fre- 
quency, the evolution of the spectral index (Fig |HJ in the 
optically-thin phase appears relatively unaffected. We take 
this as evidence that the bumps and dips in the radio light 
curve primarily reflect abrupt modulations in the CSM den- 
sity structure, rather than optical depth effects (though op- 
tical depth is tied to CSM density to some degree). We now 
consider ways in which such a structured CSM may have 
been laid down late in the life of the progenitor of SN 2001ig. 

Before doing so, we need to examine the effects of a 
change in CSM density on the velocity of the expansion, as 
well as on the radio emission. The blastwave radius increases 
with time as r a t m , where m = 1 fo r no deceleration. 
Since m = — <5/3 in the IChevalierl il982f) model, then rn = 
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0.85 for SN 200fig, implying significant deceleration in the 
surrounding CSM. The radio luminosity is related to the 
average CSM density (pcsm oc M/w) via 



L oc 



(7-7+12m)/4 



so that for SN 200fig, L oc (M/i 



(6) 



Con- 



sequently, a doubling in the CSM density will cause the radio 
emission to rise by a factor of three. 



5.1 Episodic mass-loss from a single progenitor 

The observed transition in SN 2001ig from a Type II optical 
spectrum with H lines, to a Type Ib/c spectrum without, 
argues for the ejection of a significant fraction of the red- 
giant envelope. Fig. ^indicates strong excesses of observed 
flux at t ~ 150 days, and t ~ 300 days (with a net flux ex- 
cess still at 500-600 days), hinting at a possible periodicity 
of 150 days in CSM density enhancements. If these density 
enhancements (by factors of 30% and 15% respectively) cor- 
respond to discrete shells of material expelled by the red 
supergiant, then the spacing between these shells is given 
by R s h = fexpi m . Given an initi al ejecta expansion velocity 
y e ^ of 15,000 - 20,000 km s _1 JClocchiatti fc PrietdEoOll 
IClocchiatti 120021) . an elapsed time t of 150 days, and de- 



celeration given above, then R s h — (6.3 ± 0.9) x 10 _ pc. 
Assuming the shells have been expanding at the wind veloc- 
ity w — 10 — 20 km s _1 , then the period between successive 
mass-loss episodes is T ~ 20 — 60 years. 

This is significantly longer than the timescales normally 
associated with stellar pulsations (T ~ few hundred days) 
in asymptotic giant branch (AGB) stars, but is compara- 
ble with the expected 10 2 — 10 3 yr intervals between ther- 
mal pulses (C/He she ll flashes) in 5 — 1 Mr AGB stars 
lllben fc Renzinilll983f) . Computations bv iPaczvriskil dl975h 
predict an internash period of 40 yrs when the (hydrogen- 
exhausted) core mass is nearly 1.3 Mq. This is close to the 
maximum core mass possible (depending on composition; 
iBecker fc Ibenl jl980ll 'l in stars which will finish up as a 
white dwarf, rather than undergo a supernova explosion. 
Put another way, SN 2001ig may span the gap between the 
most extreme mass-losing AGB stars, and the least-massive 
supernova progenitors. Such a massive AGB star would 
also be expected to undergo multiple convective "dredge- 
up" phases leading to modified surface ab undances, with 
N en hanced primarily at the expense of C dBecker fc Ibenl 
119791). Such changes may be apparent in the optical spectra 
jMattila et"ai]|2004h . 

There are already several good examples within our own 
Galaxy for this type of quasi-periodic mass-loss from AGB 
and post- AGB/proto- planetary nebula objects. The nearby 
carbon star IRC+10216 shows multiple nested shells spaced 
5 — 20 arcsec apart, corresponding to ejectio n timescales 
of 200 - 800 years (iMauron fc Huepins! l200Cf) . The "Egg 
Nebula" (CRL2688:ISahai et al.lll998n and IRAS 1710-3224 
jKwok. Su fc Hrivnakll998D are two bipolar proto-planetary 
nebulae which display concentric arcs with spacings on sim- 
ilar timescales. 



5.2 Binary model 

Among the other RSNe studied to date, only SN 1979C 
shows anywhere near as much systematic variation in its 
optically-thin decline phase as SN 2001ig. Indeed, for the 
first decade or so, the late-time radio light curve of SN 1979C 
could be quite well represented by a sinusoi dal modula- 
tion of the flux, with a period of 1575 days dWeiler et alJ 
1992). More recently however, these regular variations in 
SN 1979C have ceas ed, and the light curve has flattened out 
jMontes et alj EoOO). The implied modulation period of the 
mass-loss rate is T ~ 4000 years, a factor of 100 longer than 
that computed above for SN 2001ig. On the basis of evi- 
dence that the pro genitor of SN 1979C had an initial mass 
of at least 16 Mq, We iler et alJ <ll992f) argued against the 
thermal pulse scenario described in Section 15.11 as the in- 
terflash period for such a massive star and its resulting core 
would be <C 4000 years. Instead, they proposed modulation 
of the progenitor wind due to eccentric orbital motion about 
a massive binary companion as the cause of the periodicity 
in the radio light curve of SN 1979C. 

The particular binary scenario they presented had the 
red supergiant progenitor and a 10 Mq Bl dwarf orbiting 
their common barycentre, with a 4000 year period. For a 
purely circular orbit, the progenitor's orbital motion is a 
sizable fraction of the wind velocity, resulting in a spiral 
(or pinwheel-like) density structure being imprinted on the 
(otherwise uniform) mass-loss CSM in the orbital plane. 
This by itself would not lead to any periodic variation in 
the CSM density swept up by the SN shock wave. If in- 
stead the orbit was eccentric (e = 0.5 say), the acceler- 
ation of the progenitor near periastron every 4000 years 
would cause an additional pile-up of wind material which 
may then account for the observed periodicity in the radio 
emission. ISchwarz fc Pringld lll99fj) performed full hydrody- 
namical simulations of this, taking into account shocks gen- 
erated in the wind, the gravitational influence of the com- 
panion on the wind, and light-travel time effects. Their sim- 
ulations produced pronounced, but asymmetric spiral shock 
patterns, particularly when a poly tropic equation of state is 
assumed. They also demonstrated that the amplitude of the 
modulations in the radio light curve tend to decrease as our 
view of the orbital plane goes from edge-on, to face-on. This 
may partly explain why such periodic variations as seen in 
SN 1979C and in SN 2001ig are comparatively rare, since 
not only must the mass-loss rate be modulated by the right 
kind of binary orbit parameters, but we must also then be 
fortunate enough to view it from close to edge-on. 

Similar hydrodynamical simulations, but in three di- 
mensions, of detached binary systems comprising a 1.5 Mq 
AGB star and a secondary of 0. 2 5 — 2 .0 Mq were pre- 
sented by M astrodemos fc Morris! il999t) . specifically tar- 
geted at reproducing the observed structures in IRC+10216, 
CRL2688, etc., mentioned in Section 15.11 Interestingly, for 
sufficiently large binary separations (>10 AU) and wind ve- 
locities (w > 15 km s -1 ), the spiral shock structure extends 
to high latitudes, and the resulting "spiral onion shell" struc- 
ture seen in cross-section could resemble the shells seen in 
PPNe, a nd implied in SN 2001ig. In an alt ernative binary 
scenario, IHarpaz. Rappaport fc Sokei! il997f) postulate that 
it is the close passage of the secondary star effectively "chok- 
ing off" uniform mass-loss from an AGB star having an ex- 
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Figure 4. Deviations of the observed flux (on a log scale) about the best-fit model, plotted linearly as a function of time. The symbols 
represent frequencies of 22.5/18.8 GHz (circles), 15.0 GHz (crosses), 8.6/8.5 GHz (squares), 4.9/4.8 GHz (triangles), 2.4 GHz (diamonds), 
and 1.4 GHz (stars), while the solid line is a smoothed interpolation, as described in the text. 



tended atmosphere, rather than any enhancements in the 
mass-loss rate itself, which gives rise to these shells. 

Perhaps the best direct evidence for the exis- 
tence of binary-generated spiral shocks comes from 
high-resolution observations of dusty Wolf -Rayet 
(W-R) stars. iTuthill. Monnier fc Danchil (1199911 and 
iMonnier. Tuthill^^D£uichiri^99ll used the technique of 
aperture-masking interferometry on the Keck I telescope 
to image structure at better than 50 mas resolution in the 
7\-band around WR 104 and WR 98a. In both cases, they 
found pinwheel-shaped nebulae wrapping almost entirely 
around the central source out to distances of 150-300 AU. 
Their model has an OB-type companion orbiting the W-R 
star, with dust formation taking place in the wake of 
the interface region between their colliding stellar winds. 
The combination of orbital motion and wind-driven radial 
motion results in a "lawn sprinkler" effect, and the resultant 
dusty spirals. Hydrodyna mic models in 3D of the se colliding 
wind binary systems by IWalder fc FolinJ (|2003) show the 
effects of varying the orbital eccentricity. 

Remarkably, the characteristic radial scale for density 
enhancements implied by SN 2001ig's radio light curve 
(Rah = 0.0006 pc; Sect. 15. II is an almost perfect match 
to the typical scale of one full rotation of these pinwheel 
nebulae: 50-100 milliarcsec at D ~ 2 kpc -> R = 0.0005 - 
0.001 pc. Thus, SN 2001ig may represent the obliteration 



of just such a pinwheel nebula. Further support for this 
scenario comes from the flux excess (at least 60%, though 
this is a lower limit due to optical depth effects) in the first 
50 days over that expected from a simple r~ 2 CSM density 
profile, as shown in Fig. 2] Coupled with the deceleration 
mentioned previously, this would tend to favour a centrally- 
condensed additional CSM component such as a pinwheel 
nebula, rather than concentric mass-loss shells. 



5.3 SN 2001ig and the link between Type II and 
Type Ib/c SNe 

It is interesting to note that the apparent requirements to 
produce such a pinwheel nebula, i.e., a close binary sys- 
tem composed of two massive stars in which the primary 
is a W-R star, are similar to those invoked by stellar evo- 
lution theory to explain the origin of Type Ib/c SNe. The 
peculiar spectral evolution and optical light curve behaviour 
of SN 1993J and SN 1987K (Sect. [TJ has been attributed 
to t he explosion of a hydrogen-poor, helium-rich progeni- 
tor iSwartz et alJll993l) . A large fraction of the progenitor's 
original hydrogen envelope must have been shed prior to 
core collapse, either through a strong stellar wind from a sin- 
gle m assive (25-30 M@) star llHoflich. Langer fc Duschingerl 
1993); or more likely, via mass transfer from an 
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intermediate-mass (1 -15 Mp>) star in a binar y system 
jNomoto et alJ Il993j: IPodsiadlowski et all Il993t lUtrobinl 
Il994l: Van Dvk ot al. 1996). Mo dels for the evolution of mas- 
sive s t ars in close binaries (e.g., IPodsia dlowski. Joss fe Hsul 
Il992l ; IWooslev. Langer fc Weaver! Il995l) produce helium 
stars. The suggestion is that stars which have lost some or 
all of their hydrogen would be the WN class of W-R stars, 
and explode as Type lb SNe; while those which lose their 
helium layer as well would be the WC or WO classes of W- 
R stars, and explode as Type Ic S Ne llHarkness et aljll987l 
iFilfppenko. Matheson fc Holll993l) . 

At least 40% of solar neighbourhoo d W-R stars are in 
binary systems with hot companions (iMoffat et al.1 Il986t 
Ivan der Hucht et alJ Il98ftt ). and the frac tion may be even 
highe r in low-metallicity environments dPalton fc Sarazfnl 
1995), such as the outskirts of NGC 7424. The extent of mass 
transfer, and thus the end-products of the binary system, de- 
pends on the evolutionary stage of the primary at the time 
mass t ransf er commences. A s show n bv IPodsiadlowski et al] 
l|l992l) and IPols fc Nomotol l|l997T ). Case C mass transfer 
(which takes place after the core helium-burning phase) in 
systems with large eccentricity and orbital periods of a few 
years can be episodic, occurring mainly near periastron, just 
as outlined in Sect. 15.21 We propose that SN 2001ig may well 
have undergone just such a phase, without actually sharing 
a common envelope with a companion. The implied orbital 
period is given by one complete winding of the pinwheel neb- 
ula, which is simply the ratio of R B ^ ~ 0.0006 pc divided by 
the terminal wind ve locity of a WN star (~ 2000 km s ; 
lAbbott fc Contll987l) . or T ~ 100 days, consistent with this 
scenario. 

In this context, the highly-modulated radio light curves 
for SN 2001ig may represent some of the best evidence yet 
for a link between SNe of Type lib and Type Ib/c, in that 
SN 2001ig evolved optically like a Type lib, but has the radio 
characteristics that should be expected for a Type Ib/c SN 
originating in a W-R + OB binary system viewed nearly 
edge-on. A testable prediction from this scenario is that 
the companion star (which by virtue of mass accretion may 
be even brighter than the W-R progenitor of SN 2001ig) 
should eventually become visible against the fading opti- 
cal remnant, a s has been postulated but not yet observed 
for S N 1993J iPodsiadlowski et al]ll993 iFilippenko et al] 
Il993ft . 

Besides its apparent association with the gamma-ray 
burst GRB 980425, SN 1998bw is also notable for show- 
ing bumps and dips in its radi o light curves not di ssimi- 
lar to those seen in SN 200 lig iWeiler. Panaeia fc Montesl 
12001]) . These deviations, while exhibiting no clear period- 
icity in the first 1000 days, were less conspicuous at low 
frequencies, just as in SN 2001ig. This kind of behaviour 
is likely due to the CSM still being optically thick to low 
frequencies at relatively late times, and therefore only emis- 
sion originating in the near-side of the CSM can be seen; 
whereas at high frequencies, emission from the entire CSM 
is visible. SN 1998bw had the spectral characteristics pri- 
marily of a Type Ic e vent, with a probable W-R progenitor 
ilwamoto et al] I1998T) . Unfortunately, barely a half-dozen 
Type Ib/c SNe have been detected or studied in the radio, 
so it is too early to tell whether such modulated radio light 
curves may be the clue to a common massive-binary origin 
for Type lib and Type Ib/c SNe, and possibly some GRBs. 



6 CONCLUSIONS 

By compiling one of the most complete multi-frequency ra- 
dio datasets ever collected for any supernova, we have been 
fortunate to witness regular modulations in the radio light 
curves of the Type lib Supernova 2001ig. The time taken 
to reach peak luminosity at 5 GHz, the rate of decline since 
then, and the derived mass-loss rates prior to explosion are 
all intermediate between those of Type Ib/c and "normal" 
Type II SNe. We find the light curve modulations to recur 
on a timescale of ~ 150 days, and have shown them to be 
true density modulations in the CSM, and not optical depth 
effects. Allowing for the deceleration of the ejecta, these den- 
sity enhancements are spaced 0.6 milli-parsecs (or 130 AU) 
apart. 

While we cannot totally exclude the possibility that 
these density enhancements represent mass-loss shells from 
the thermal-pulsing phase of a single AGB star progenitor, 
we find the weight of evidence supports a stellar wind, mod- 
ulated by motion in an eccentric binary system, as their 
source. As had been suggested previously for SN 1979C, the 
combination of a massive binary companion causing a pile- 
up of mass-loss during periastron, and a favourable viewing 
angle, can result in just the kind of periodic density varia- 
tions observed in SN 2001ig. Recent near-IR interferometric 
observations of the anticipated "pinwheel" dust nebulae in 
systems comprising a Wolf-Rayet star and a hot massive 
companion lend weight to this scenario, especially as the 
observed size scales match those required for SN 2001ig. Fi- 
nally, as optical spectroscopy has recently been leading to 
the conclusion that Type Ib/c SNe are the product of a core- 
collapse event in the W-R component of such systems (and 
Type lib SNe being those caught early enough to reveal the 
last traces of their lost hydrogen envelope) , we believe these 
radio observations of SN 200 lig provide the "missing link" 
between the pinwheel nebulae observed in Galactic W-R bi- 
nary systems, and their eventual fate in Type lib or Type 
Ib/c superriovae. 
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